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I.  INTRODUCTION 


OC  various  methods  proposed  for  the  removal  of  a  chemical  contaminant 
from  s.irt'aces  of  vehicles  or  other  equipment  of  Army's  Interest,  utilization 
of  liquid  jet  spray  is  the  most  practical  and  effective  one  at  the  current 
levei  >f  technological  development .  It  uses  the  force  produced  by  the 
turnin''  of  jet  stream  at  the  impingement  to  displace  the  contaminant.  The 
surfaces  then  can  be  decontaminated  by  moving  the  jet  stream  toward  the 
co  ita  ni aunt . 

The  chemical  contaminant  is  deployed  on  the  surfaces  in  the  form  of 
droplets  with  .a  number  density  of  approximately  420  droplets/m^  .  Each 
droplet  has  a  diameter  of  2  to  4  ram.  The  density  oi  the  contaminant  is  very 
close  to  that  of  water,  however,  its  viscosity  can  be  from  10  to  1000  times 
higher  than  water  viscosity,  depending  on  the  enviromental  condition. 

A  high-performance  jet  spray  for  decontamination  should  possess  the 
following  two  features:  a  high  cleaning  speed  and  an  efficient  use  of  jet 
fluid.  These  features  are  particularly  important  when  the  jet  system  is 
operated  in  the  field  where  it  is  often  required  to  decontaminate  a  surface 
area  in  the  shortest  period  of  time  and  with  the  least  consumption  of  jet 
fluid.  In  designing  such  a  system,  knowledge  of  fundamental  characteristics 
of  the  flow  interaction,  such  as  the  evolution  of  the  contaminant  droplet 
and  the  effects  of  various  flow  parameters  on  the  flow,  is  vital.  This 
information  can  be  sought  via  experiments,  but  that  would  require  very 
sophisticated  instrumentation  and  would  be  very  costly.  As  an  alternative, 
one  may  use  computer  simulations  based  on  appropriate  flow  models.  This 
method  can  have  a  much  greater  flexibility  than  experiments  for  examining 
areas  of  importance  in  the  flow  field  and,  therefore,  can  provide  better 
insights  into  the  flow  phenomena. 

The  interaction  flow  involves  two  prime  fluids,  the  jet  fluid  (water) 
and  the  contaminant,  separated  by  Interfaces.  They  also  may  have  free 
sufaces  with  the  ambient  (air)  .  The  problem  is  further  complicated  by  the 
fact  that  the  flow  is  three-dimensional  and  is  highly  transient .  Msst  of 
the  investigations  undertaken  in  the  past  in  jet  impingements  relate  to  the 
VTOL  program  (vertical  takeoff  and  landing  aircraft)  or  rockect  exhaust 
flows,  and  they  are  concerned  with  jet  impingements  on  a  solid  plane. 

‘ G .  T.  Thy  lor,  "Oblique  Impact  of  a  Jet  on  a  Plane  Surface ,  "  Phil .  Trane ,  R. 
Soc.  A,  260,  1966,  pp.  96-100.  ‘ 

2</.  H.  Michell,  Phil.  Trans.  A,  181,  1890,  pp.  S89-431. 
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A.  Rubell,  "Computations  of  Jet  Impingement  on  a  Plat  Surface,  "  AIAA  J., 
18,  No.  2,  February  1980,  pp.  168-175.  ~~ 

lA.  Rubell,  "Computations  of  the  Oblique  Impingement  of  Round  Jets  upon  a 
Plane  Vail,  "  AIAA  J.,  19,  No.  7,  July  1981,  pp.  863-871. 


Steady  two-dimensional  jet  impingements  on  a  liquid  surface  were  considered 
by  Hunt"*  and  Vanden-Broeck.^  They  characterized  the  wave-like  hydrodynamic 
Instability  that  occurs  at  the  interface  of  a  gas  jet  and  a  liquid  surface 
by  using  simplified  theories,  but  made  no  predictions  of  the  movement  of  the 
liquid  and  the  pressure  distribution  in  the  flow  field,  which  are  of 
Interest  to  us. 

In  this  study  we  focus  our  attention  on  the  interaction  of  a  water  jet 
with  a  single  contaminant  droplet.  We  simplify  the  problem  by  treating  it 
as  a  two-dimensional  flow.  Two  models,  called  one-fluid  flow  and  two-fluid 
flow,  were  developed  to  treat  two  pre-impingement  configurations:  present 
and  absent  water  coverage  over  the  droplet.  The  flow  field  is  governed  by 
the  unsteady  Navier-Stokes  equations  which  were  solved  numerically  via 
finite  difference  schemes  using  the  SOIA-VOF  code.^  We  implemented  a 
viscosity  relation  in  the  computer  code  to  adapt  it  to  the  present  case  with 
two  fluids  having  distinct  viscosities. 

The  purpose  of  the  present  study  is  to  characterize  the  Jet-contaminant 
interaction  and  to  provide  useful  data  for  design  of  the  jet  spray.  Results 
presented  in  this  report  consist  of  computer-generated  flow  patterns  which 
show  the  flow  development,  effects  of  various  parameters  on  the  flow,  and 
Impact  pressures  on  the  impingement  wall.  Also  presented  are  flow 
developments  of  jet-contaminant  interactions  in  cavities  and  at  corners  of 
two  perpendicular  walls. 

II.  FLOW  MODELS 

Figure  1  depicts  two  pre-impingement  flow  configurations  which  can 
occur  in  decontamination  processes.  In  the  first  configuration,  shown  in 
Figure  la,  a  water  jet  is  directed  at  a  contaminant  droplet  with  a  wacer 
coverage  over  it  on  a  plane  wall.  The  second  configuration,  shown  in  Figure 
lb,  differs  from  the  first  in  that  there  is  no  initial  water  coverage  over 
the  droplet.  In  order  to  characterize  the  flows  developing  from  these  two 
configurations,  we  have  developed  a  two-fluid  model  and  a  one-fluid  model. 
Both  models  describe  a  two-dimensional  viscous  flow  and  can  be  handled  by 


5</.  N.  Hunt ,  'Mate  Formation  in  Explosive  Welding,  "  Philosophical  Magazine, 
8th  Seriee,  Vol.  17,  p.  669-680,  1967 . 

/? 

J.  M.  Vanden-Broeck,  "Deformation  of  a  Liquid  Surface  by  an  Impingement  Gae 
Jet,  "  SIAM  J.  Appl.  Math.,  41,  No.  2,  October  1981,  pp.  306-309. 

?B.  D.  Nicholas,  C.  W.  Hurt,  and  R.  S.  Hotchkiss,  "SOLA-VOF:  A  Solution 
Algorithm  for  Transient  Fluid  Flow  with  Multiple  Free  Boundaries,  "  Los 
Alamos  Scientific  Laboratory  Report  No.  LA-8355,  1980. 


^  PLANE  WALL 

a.  Droplet  with  Water  b.  Droplet  without  Water 

Layer  Coverage  (Two-  Coverage  (One-Fluid 

Fluid  Model)  Model) 

Figure  1.  Pre-impingement  Flow  Configurations 

the  SOLA-VOF  computer  code  which  can  handle  flow  problems  involving  two 
immiscible  fluids  separated  by  interfaces  in  a  region  without  voids  or  flow 
problems  involving  one  fluid  with  voids  (air). 

2.1  TVo— Fluid  Model 

This  is  a  model  of  a  channel-type  flow  bounded  by  the  dashed  line 
indicated  in  Figure  la,  covering  the  major  part  of  the  flow  region.  As 
shown  in  Figure  2,  the  channel  contains  two  fluids  (water  and  contaminant) 
separated  by  an  Interface.  The  upper  wall  of  the  channel  coincides  with  the 
upper  boundary  of  the  water  layer  so  as  to  eliminate  consideration  of  its 
free  surface  with  the  ambient.  An  outflow  boundary  condition  is  specified 
at  this  wall  and  at  the  ends  of  the  channel,  allowing  the  fluids  to  flow  out 
the  region.  The  contaminant  which  occupies  the  shaded  rectangular  region  is 
assumed  to  perfectly  wet  the  lower  wall  of  the  channel.  To  account  for 
viscous  effects,  a  no-slip  condiation  is  used  at  this  wall.  Finally,  a 
steady  and  uniform  jet  stream  at  an  angle  of  incidence,  ©  ,  is  prescribed 
along  a  segment  of  the  upper  wall. 

2.2  One-Fluid  Model 

Results  of  computations  with  the  two-fluid  model,  given  in  Section  5-B 
of  this  report,  show  that  for  a  close-in  impingement  the  movement  of  the 
contaminant  droplet  at  its  early  stage  is  insensitive  to  the  the  viscosity 
of  jet  fluid.  In  addition,  the  density  of  the  contaminant,  1070  kg/nr*,  is 
very  close  to  that  of  water,  1000  kg/nf*.  Therefore,  we  can  simplify  the 
problem  by  setting  the  physical  properties  of  the  jet  fluid  (water)  equal  to 
that  of  the  contaminant  and  treat  the  jet-contaminant  interaction  as  a  one- 
fluid  flow  problem.  Figure  3  depicts  the  flow  region  for  computation, 
bounded  by  the  dashed  line  marked  in  Figure  lb.  In  this  model,  only  the 
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free  surface  with  the  ambient  needs  to  be  traced,  but  not  the  water- 
contaminant  Interface  since  both  liquids  are  assumed  to  be  Identical.  The 
model  Is  suitable  for  the  treatment  of  the  flow  configuration  of  Figure 
lb.  The  two-fluid  model  cannot  be  used  for  this  configuration  because  It 
would  require  the  tracing  of  both  the  free  surface  and  the  interface,  and 
thus  would  be  beyond  the  capability  of  the  SOLA-VOF  code. 


III.  FLOW  EQUATIONS,  KTHOD  Of  SOLUTION,  AND  MSH  SETUP 


3.1  Flow  Equations  and  Ifethod  of  Solution 


The  governing  equations  of  the  above  two  model  flows  are  the  continuity 
equation 
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and  the  momentum  equations 
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where  t  is  time  u  and  v  are  the  x-component  (along  the  channel)  and  the  y- 
component  (normal  to  the  channel)  of  the  flow  velocity,  respectively.  The 
density  p,  the  sound  speed  c,  and  the  kinematic  viscosity  v  are  assumed  to 
be  constant.  The  Reynolds  numbers  based  on  the  jet  thicknesses  and  the  jet 
velocities  used  in  this  study  are  between  20  and  2000.  Within  this  range, 
Eqs .  (2)  and  (3)  are  felt  to  be  appropriate  without  considering  turbulence 
effects . 

In  order  to  track  the  water-contaminant  interface  in  two-fluid  flows  or 
the  free  surface  in  one-fluid  flows,  the  SOLA-VOF  code  uses  a  "fractional 
volume  of  fluid"  function  F.  The  function  satisfies  the  following  relation 


which  states  that  F  moves  with  the  fluid.  In  two-fluid  flows,  the  value  of 
F  in  a  computational  cell  is  equal  to  the  fractional  volume  of  the  cell  that 
is  occupied  by  the  first  fluid  (say,  contaminant) .  A  unit  value  of  F 
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corresponds  Co  a  cell  full  of  the  first  fluid,  whereas  a  zero  value 
Indicates  that  the  cell  contains  only  the  second  fluid  (water  from  the  jet 
and  the  water  layer  over  the  contaminant) .  Cells  with  F  values  between  zero 
and  one  contain  an  Interface,  as  Illustrated  In  Figure  4.  In  one-fluid 
flows,  the  second  fluid  has  the  meaning  of  a  void  and  F  traces  the  free 
surface . 


Figure  4 .  Free  Surface  (or  Interface) 

Across  a  Mesh  Cell 

The  SOLA-VOF  Code  solves  the  momentum  Eqs .  (2)  and  (3)  by  using  an 
explicit  finite  difference  scheme  and  the  continuity  Eq.  (1)  by  using  an 
Implicit  finite  difference  method.  The  solution  of  Eq .  (4)  is  obtained  by 
using  a  Donor-Acceptor  flux  approximation.  Details  of  the  solution  method 
are  given  In  Reference  7 .  In  order  to  better  observe  the  evolution 
(location  and  shape)  of  the  droplet,  marker  particles  were  embedded  in  the 
region  originally  covered  by  the  droplet.  These  particles  move  with  the 
fluid,  but  do  not  affect  the  flow. 

The  SOLA-VOF  code  requires  that  the  viscosities  of  the  two  fluids  In  a 
two-fluid  flow  case  are  equal .  To  adapt  the  code  to  the  present  problem 
involving  two  fluids  with  distinct  viscosities,  we  used  the  viscosity 
relation 


v  *  v  F  +  (1-F)  v 

c  w 


(5) 


16 


where  v  and  v  are  the  kinematic  viscosities  of  contaminant  and  water, 
respectively,  "and  v  is  the  average  kinematic  viscosity  of  the  two  fluid 
mixture  in  the  cell.  Similarly,  the  density  in  a  cell  was  approximated  by 


p  -  Pc  F  +  (1-F)  pw  (6) 

By  Eqs .  (5)  and  (6)  the  values  of  v  and  p  in  a  cell  are  averages  weighed  by 
the  value  of  F.  For  instance,  in  a  cell  full  of  contaminant  the  value  of  F 
is  one  and  thus  v  ■  vc  and  p  ■  pc  . 
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For  the  numerical  computations  an  Eulerian  mesh  was  used,  consisting 
of  rectangular  cells  of  variable  sizes  with  finer  zoning  near  the 
impingement  wall .  The  use  of  a  variable  mesh  can  substantially  save 
computer  costs  without  sacrificing  the  accuracy  of  computation  in  the 
regions  of  main  Interest. 


3.2.1  TWo— Fluid  Flow .  Figure  2  shows  part  of  the  100x12  mesh  of  the 
flow  region  to  be  analysed.  The  mesh  size  in  the  x-dlrectlon  is 
uniform,  Ax**0 .2  mm,  while  the  mesh  size  in  the  y-directlon  is  gradually 
reduced  from  Ay=0 .1  mm  at  the  top  to  Ay**0 .025  mm  (approximately  4Z  of  the 
droplet  height)  at  the  bottom  of  the  flow  channel.  In  order  to  test  the 
appropriateness  of  the  mesh,  computations  with  slightly  finer  as  well  as 
slightly  coarser  meshes  were  performed .  The  only  visible  variation  noticed 
was  the  location  of  the  tip  of  the  upstream  edge  of  the  droplet  immediately 
above  the  wall.  This  variation  is  to  be  expected  because,  as  illustrated  in 
Figure  5,  the  mesh  size  near  the  wall  directly  affects  the  calculated 
magnitude  of  the  displacement  S  of  the  droplet  upstream  edge .  This  effect 
will  be  discussed  in  Subsection  A-l  of  Chapter  5. 


Figure  5  .  Relation  Between  Mesh  Size  and 

Displacement  S  of  Droplet  Upstream  Edge 


3.2.2  One-Fluid  Flow.  Figure  3  depicts  the  computational  domain  for 
the  one-fluid  flow  analysis.  Initially,  a  uniform  velocity  is  specified  in 
all  mesh  cells  inside  the  jet  region  and  the  velocity  is  zero  elsewhere. 
For  subsequent  time  levels,  only  the  mesh  cells  along  the  upper  boundary  of 
the  jet  region  supply  a  steady  and  uniform  jet  stream. 

The  lower  part  of  Figure  3  shows  a  portion  of  the  50x18  mesh  for 
numerical  computation.  The  mesh  size  was  chosen  for  the  purpose  to  generate 
patterns  using  the  one-fluid  model.  Determination  of  effects  of  flow 
parameters,  which  requires  a  higher  accuracy,  was  done  by  using  the  two- 
fluid  model. 


IV.  TEST  RDN  OF  COMPUTER  CODE  AND  CONSIDERATIONS 
OF  SURFACE  TENSION  AND  CONTACT  ANGIE 

4.1  Test  Run  of  the  SOLA-VOF  Code 

As  a  test  run  of  the  computer  code,  we  calculated  the  pressure 
distribution  on  a  plane  wall  impinged  by  to  a  two-dimensional, 
incompressible,  inviscid,  normal  jet  impingement.  The  result  was  compared 
with  the  steady-state  pressure  calculated  from  the  following  relations  given 
in  Tayor's  paper: * 

P(Xl)  -  (1/2)  pvj2(l-u^)  (7) 

Uj  =  (-1  -  q  cos  0  )  +  /  ( 1  -  qZ)  sin  0  /  (q  —  cos  G  )  (8) 

Xj*  (1/2) 1(1+  cos  0)  ln( 1  +  q)  -  (1  -  cos  0)  ln( 1  -  q)] 

+  sin  0  sin-*  q  (9) 

where  u  ■  U/Vj  (U  being  the  fluid  velocity  on  the  wall),  Xj  **  tt  x/(thickness 
of  jet),  x  is  the  distance  along  the  wall  from  the  stagnation  point,  q  is  an 
auxiliary  variable  to  be  eliminated  between  Eq.  (8)  and  Eq.  (9),  and  6  is 
the  angle  of  jet  incidence. 

Figure  6  is  the  setup  of  flow  region  for  the  test  problem.  There  Is  a 
fluid  column  initially  standing  on  the  plane  wall  with  a  uniform  flow 
velocity  100  m/s  specified  in  the  fluid  region.  For  subsequent  time  levels, 
the  uniform  flow  velocity  remains  only  along  the  upper  boundary  of  the 
column,  while  the  lower  part  of  the  column  deforms  continuously  until  a 
steady  flow  is  established.  Plots  in  Figure  7  show  the  flow  development 
generated  by  the  computer.  The  instantaneous  pressure  profiles  on  the 
impingement  wall  at  times  up  to  0.4  millisecond  (ms)  are  presented  in  Figure 


8.  AC  tine  equal  Co  0.4  ns,  Che  pressure  approaches  Che  steady-state  valve 
and  Che  pressure  profile  is  very  close  Co  the  steady-state  pressure  profile 
calculated  from  Eq.  (7). 

4.2  Surface  Tension 

The  surface  tension  at  the  water-contaminant  interface  has  the 
magnitude  of  1-2  dynes/mm. ®  The  effect  of  the  surface  tension  was  tested  by 
computing  the  displacement  S  of  the  droplet  upstream  edge  versus  time  for  a 
sample  flow.  The  resulting  displacements  and  droplet  profiles  are  given  in 
Table  1  and  Figure  9,  respectively.  There  is  virtually  no  difference 
between  results  with  and  without  consideration  of  surface  tension. 
Therefore,  in  order  to  save  computer  time,  we  neglected  the  surface  tension 
in  subsequent  computations. 
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Figure  6.  Setup  for  a  Normal  Jet  Impingement 
on  a  Plane  Wall 

Table  1.  Displacement  S  (mm)  of  Droplet  Upstream  Edge  vs.  Time 

after  Commencement  of  Jet  Flow  for  Cases  with  and  without 
Considerations  of  Surface  Tension  and  for  Different 
Contract  Angles 


TIME 

SURFACE  TENSION 

- a _ _ 

I  zaynes/mm  1 

IHI5I 

CONT.  ANG.*90 

CONT.  ANG.  *45° 

CONT.  ANG.  *90u 

0.04 

0.195 

0.195 

0. 195 

0.08 

0.428 

0.427 

0.428 

0. 12 

0. 673 

0.673 

0.673 

0.16 

0.910 

0.910 

0.910 

0.20 

1.149 

1.147 

1.148 

0.24 

1. 369 

1.369 

1.370 

0.28 

1.583 

1.582 

1.583 

V.-  5  m/s,  v-  9. 8  mm*/s,  8  -  56. 28° 
J  c 


J.  Matta,  U.S.  Army  Chemical  Research  and  Development  Center  ( CRDC), 
private  communication,  January  1983. 


4  .3  Contact 


Because  the  contact  angle  of  the  contaminant  with  the  wall  surface  has 
not  been  well  documented,  It  was  assumed  to  be  90°  In  this  study.  To  test 
the  significance  of  the  angle  for  the  results,  we  picked  another  angle  of 
45°  for  a  sample  calculation.  Table  1  lists  the  results  in  comparison 
to  the  90°  case.  The  sample  calculation  shows  that  the  contact  angle  does 
not  Influence  the  results . 


t  ■  0. 1  ms 
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t  ■  0. 4  ms 


0.40 

CALCULATED  FROM 
EQ.  (7) 


V.  Cdf*  OTATIONAL  RESULTS 

We  present  results  for  the  following  ranges  of  input  data: 

9  *  angle  of  jet  Incidence  ■  0°  -  90° 

Dj  ■  jet  size  (thickness  of  jet  in  two-dimensional  model) 

■  1 .83  -  3 .66  mm 

Vj  =  jet  velocity  ■  5  -  12.5  m/s,  uniform  and  steady 

Pw  *  water  density  -  1000  kg/nr* 

=  contaminant  density  ■  1070  kg/m^ 

=  kinematic  viscosity  of  water  ■  0.98  mm^/s 
vc  -  kinematic  viscosity  of  contaminant  ■  9.8  -  980  mm^/s 

Droplet  dimensions:  3  mm  x  0.6  mm  (length  by  height) 


The  jet  velocities  chosen  produce  stagnation  pressures  of  12.4  -  75.8  kPa 
(1.8  -  11  psi)  ,  which  are  practical  for  decontaminations.  The  dimensions  of 
the  contaminant  droplet  represent  the  average  droplet  size  deployed  on  a 
plane  wall.  Our  computational  results  indicate  that  for  such  a  droplet  size 
the  jet  thicknesses  given  above  can  perform  the  decontamination  effectively 
and  efficiently. 

The  principal  results  presented  include: 

a)  flow  patterns 

b)  contaminant  viscosity  effects 

c)  optimum  angle  of  jet  incidence  at  which  the  jet  performs  best 

d)  effects  of  jet  velocity  and  jet  size  on  jet  performance 

e)  improvements  of  decontaminating  speed  and  efficiency  of  jet  fluid  usage 

f)  pressure  distribution  on  the  impingement  wall 


Some  of  the  results  were  only  obtained  by  using  the  two-fluid  model  which 
has  a  higher  accuracy.  These  results  are  presented  in  the  next  section. 
Results  obtained  by  the  one-fluid  model  will  be  presented  in  Section  5-B. 
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5.1  TWo-Fluld  Model 


5.1.1  Flow  Patterns  and  Viscosity  Effects.  Figure  10  presents  two  sets 
of  flow  development  for  times  up  to  0.2  ms  after  the  start  of  the  jet 
impingement.  They  correspond  to  the  two  contaminant  viscosities  v c  *  9.8 
mnr/s  (10  v  )  and  =980  mm^/s  (1000  v  w) ,  respectively.  The  plots  have 
been  magnified  by  a  factor  of  3  in  the  vertical  direction  in  order  to 
provide  a  clear  flow  visualization  near  the  bottom  wall.  The  arrows  in  the 
flow  region  represent  the  local  fluid  velocities.  The  figures  show  that  the 
contaminant  viscosity  has  a  pronounced  effect  on  the  displacement  S  of  the 
droplet  upstream  edge  as  well  as  on  the  droplet  profile.  In  the  low 
viscosity  case  the  jet  stream  toward  the  droplet  is  nearly  parallel  to  the 
bottom  wall,  whereas  in  the  high  viscosity  case  the  jet  stream  is  lifted  off 
the  wall.  It  is  interesting  to  note  that  as  the  time  progresses,  part  of 
the  fluid  near  the  downstream  edge  of  the  droplet,  especially  in  the  high 
viscosity  case,  moves  toward  the  bottom  wall  behind  the  droplet. 
Apparently,  a  low  pressure  region  is  created  behind  the  lower  portion  of  the 
droplet,  similar  to  the  flow  phenomena  that  appears  behind  a  stationary 
obstacle.  Figure  11  is  a  map  of  the  marker  particles  embedded  in  the 
droplet  region. 

The  displacement  S  indicated  in  Figure  11  is  defined  as  the  distance 
from  the  initial  location  of  the  upstream  edge  of  the  droplet  (marked  by  the 
dashed  line)  to  the  nearest  marker  particle  on  the  first  row  above  the 
inpingement  surface.  We  embedded  24  rows  of  marker  particles  uniformly 
across  the  droplet  in  the  vertical  direction,  and  there  were  30  marker 
particles  in  each  row.  As  illustrated  in  Figure  12,  the  nearest  marker 
particle  from  the  dashed  line  is  not  always  the  one  which  originally  resided 
at  the  very  left  end  of  the  first  row.  This  is  because  the  jet  stream  is 
lifted  off  the  surface  and,  as  a  consequence,  the  marker  particles  in  the 
front  are  carried  upward.  The  location  of  each  marker  particle  following 
the  jet  impingement  can  be  traced  and  printed  out  by  the  computer.  The 
displacement  S,  as  defined,  is  dependent  on  the  number  of  rows  embedded  and 
on  the  fineness  of  the  computational  grid.  A  larger  number  of  rows  and  a 
finer  mesh  size  would  result  in  a  smaller  displacement  S  and  theoretically 
S  ♦  0  as  the  marker  particle  approaches  the  wall.  In  the  present  study,  the 
finest  cell  size  In  the  vertical  direction  used  near  the  surface  was  0.025 
mm,  which  is  approximately  4%  of  the  droplet  height.  The  initial  distance 
between  the  first  row  of  marker  particles  and  the  surface  was  also  0.025 
mm.  Hence  our  definition  of  S  is  essentially  the  displacement  of 
contaminant  particles  that  are  initially  more  than  0.025  mm  off  the  wall. 
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5.1.2  C^tlarn  Angle  of  Jet  Incidence.  The  angle  of  jet  Incidence 
plays  a  key  role  In  displacing  the  contaminant.  A  small  angle  causes  a 
large  fraction  of  jet  stream  to  move  toward  the  droplet,  but  a  small 
impingement  force  on  the  contaminant  near  the  wall.  A  large  angle,  on  the 
other  hand,  has  reverse  effects.  Therefore,  there  is  an  optinum  angle  at 
which  a  jet  performs  best  in  removing  the  contaminant. 

Prior  to  proceeding  to  optimization  of  the  angle  of  incidence,  we  have 
to  note  that  the  performance  of  the  jet  spray  should  be  evaluated  from  two 
aspects,  namely,  the  cleaning  speed  and  the  efficiency  of  jet  fluid  usage. 
The  efficiency  of  jet  fluid  usage  may  be  important  in  the  areas,  especially 
in  the  field,  where  the  supply  of  jet  fluid  can  be  limited. 


In  order  to  optimize  the  angle  of  jet  incidence,  we  used  several  values 
between  19.8°  and  90°  to  compute  the  displacement  of  the  droplet  as  a 
function  of  time  following  the  impingement.  At  each  angle  the  jet  was 
located  such  that  the  impingement  takes  the  shortest  time  to  displace  the 
upstream  edge  of  the  droplet  a  prescribed  distance,  one  third  of  the 
original  droplet  length.  A  larger  prescribed  distance  likely  would  not 
change  the  results  because  the  movement  of  the  droplet  at  the  very  early 
stage  suffices  to  define  the  performance  of  the  jet.  The  angle  was 
optimized  by  maximizing  the  cleaning  speed  and  minimizing  the  expenditure  of 
jet  fluid.  The  particular  values  of  the  angle  for  the  computations  were 
chosen  to  be  conveniently  adapted  to  the  computational  mesh. 

Cleaning  Speed 


The  cleaning  speed,  which  is  the  rate  of  decontamination,  is  defined  as 
an  average  velocity  of  the  movement  of  the  droplet  upstream  edge  along  the 
wall.  The  average  is  taken  over  the  time  t  needed  to  move  the  edge  by  S  mm, 
that  is  S  **  S/t.  Figures  13  and  14  show  the  displacement  S  of  the  droplet 
upstream  edge  versus  time  for  the  jet  velocities  V.  -  10  m/s  and  V  ■  5  m/s, 
respectively.  The  dashed  lines  represent  the  results  for  the  case  that  the 
contaminant  viscosity  is  equal  to  10  times  water  viscosity  (i.e.,  *  10  vw  ) 
and  the  solid  lines  for  V£  *  1000  v  .  Figures  15  and  16  present  the 
corresponding  mean  velocities  and  show  that  S  at  a  given  time  t  increases 
with  the  angle  of  incidence,  9,  until  9  reaches  56.28°.  Beyond  that  angle, 
S  decreases  except  in  the  case  of  low  jet  velocity  and  extremely  large 
viscosity  shown  in  Figure  16.  When  the  water  layer  which  initially  covers 
the  droplet  shown  in  Figure  la  is  increased  from  0.2  mm,  which  is  used  for 
Figures  13  and  14,  to  0.6  mm,  the  result  shown  in  Figure  17  also  indicates 
that  9  *  56.28°  results  in  a  better  jet  performance  than  any  other  angle. 
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When  the  contaminant  viscosity  is  extremely  large  and  the  jet 
velocity  is  low,  the  normal  force  is  the  principal  mechanism  to  break  up  and 
displace  the  droplet.  Because  the  normal  force  is  produced  by  the  turning 
of  flow  direction  and  is  proportional  to  8,  S  increases  monotonlcally 
with  8  to  90°  in  such  cases,  as  shown  by  the  solid  curves  in  Figure  14. 

Efficiency  of  Jet  Fluid  Osage 

In  Figure  18,  "A"  denotes  the  surface  area  which  has  been  cleaned 
(decontaminated)  and  "V*"  the  volume  of  jet  fluid  consumed  for  cleaning  the 
area.  For  a  unit  width  of  jet  in  the  present  two-dimensional  model,  we 
define  that  A  ■  S  **  S  t  m^  and  V*  =  Vj  Djt  m^.  Then  A/V*,  called  the 
nominal  area  cleaned  per  unit  volume  of  jet  fluid  consumed,  indicates  how 
efficiently  the  jet  fluid  is  used.  The  curves  show  that  a  jet  at  an  angle 
of  incidence  between  45°  and  60°  will  achieve  the  most  efficient  use  of  jet 
fluid.  They  also  show  that  when  8  is  reduced  below  40°,  the  efficiency 
falls  sharply. 


5.1.3  Effects  of  Jet  Velocity  on  Jet  Performance.  We  expect  that  the 
jet  velocity  affects  both  the  cleaning  speed  and  the  efficiency  of  jet  fluid 
usage.  A  larger  jet  velocity  produces  a  greater  normal  force  on  the 
contaminant  and,  accordingly,  more  cleaning  power,  as  shown  in  Figure  19. 
This  result  can  also  be  recognized  when  we  compare  Figure  16  with  Figure  15, 
in  which  the  cleaning  speed  has  been  Increased  substantially  as  the  jet 
velocity  Vj  is  raised  from  5  m/s  to  10  m/s.  A  larger  jet  velocity,  however, 
consumes  more  jet  fluid  in  a  given  period  of  time  and  its  Impact  on  the 
efficiency  of  jet  fluid  usage  has  to  be  examined. 

Figure  20  reveals  that  the  efficiency  monotonlcally  Increases  with  the 
jet  velocity  in  the  range  studied. 


5.1.4  Effects  of  Jet  Size  on  Jet  Performance.  In  the  two-dimensional 
model,  the  jet  size  Dj  is  simply  the  thickness  of  the  jet.  If  Dj  is  too 
small,  say,  less  than  the  height  (0.6  mm)  of  the  droplet  on  a  plane  wall, 
the  jet  stream  towards  the  droplet  can  not  cover  the  entire  leading  surface 
of  the  droplet  and  it  becomes  a  jet  flow  submerged  in  the  droplet  as 
depicted  in  Figure  21.  The  cleaning  power  will  certainly  be  low.  Although 
no  extensive  confutations  have  been  performed  in  this  study  to  determine  the 
minimum  size  required  for  D..,  indications  based  on  the  results  in  various 
flow  cases  show  that  a  value  larger  than  one  half  the  droplet  length 
(diameter  in  three-dimensional  cases)  is  adequate. 

Figure  22  reveals  that  the  displacement  S  at  a  given  time  and 
consequently  the  cleaning  speed  increase  with  jet  size  Dj.  However,  the 
efficiency  of  jet  fluid  usage  indicated  in  Figure  20  decreases  gradually  as 
Dj  is  increased  in  both  low  and  high  contaminant  viscosity  cases. 


1*3* 


U 

5 

S50.75 


soa  h  It3 


B  -  56. 28 
Dj  •  1.83  mm 


- U-YOU 


0  0.(5 


0.1  0.15 

TIME,  t  ms 


Figure  19.  Displacement  of  Droplet  Upstream  Edge,  S,  vs.  Time  after 
Commencement  of  Jet  Flow,  t,  for  Various  Jet  Velocities,  V 


5.1.5  An  Kffectlve  and  Efficient  Means  to  Improve  Jet  Perf organ 

We  have  shown  in  the  previous  two  subsections  that  an  increase  of  either 
jet  velocity  or  jet  size  can  improve  the  cleaning  speed.  Now  we  investigate 
which  is  more  advantageous  to  improve  the  cleaning  speed  for  a  given  jet 
flow  rate  (6  -  Vj  Dj)  by  increasing  the  jet  velocity  or  by  increasing  the 
jet  size.  To  evaluate  the  computed  results  we  first  compare  curve  "a"  with 

curve  "a*”  in  Figure  23  or  Figure  2  4.  As  listed  at  the  right  corners  of^the 

figures,  both  curves  correspond  to  the  same  flow  rate  Q  ■  0.0183  nr/s, 
however,  curve  "a"  is  for  a  jet  velocity  twice  that  of  curve  "a*”  and  a  jet 
size  half  that  of  curve  "a*".  A  comparison  between  them  reveals  that  a 
higher  jet  velocity  produces  a  higher  cleaning  speed.  This  is  also  true 
when  we  conpare  curve  "a”  with  "a*",  curve  "b"  with  curve  "b*”,  or  curve  "5" 
with  curve  "b*”,  all  for  different  combinations  of  Vj  and  Dj.  In  the 
previous  two  subsections  and  Figure  20,  we  also  have  seen  that  the 

efficiency  of  jet  fluid  usage  can  be  raised  by  increasing  the  jet  velocity. 
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Figure  21.  Flow  Developing  from  an  Insufficient  Jet  Size 
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but  not  by  enlarging  the  jet  size.  Therefore,  It  is  more  advantageous  to 
increase  the  jet  velocity  than  the  jet  size.  An  inportant  consequence  for 
practical  use  is  that  a  jet  spray  composed  of  a  number  of  small  high-speed 
jets  is  more  effective  and  more  efficient  than  a  single  large  low-speed  jet 
with  the  same  flow  rate. 


5.1.6  Pressure  Distribution  on  Impingement  Wall.  Figure  25  shows  some 
instantaneous  pressure  distributions  on  the  impingement  wall.  The  peak 
pressure  can  be  more  than  11  times  the  corresponding  steady-state  stagnation 
pressure  of  the  jet,  (l/2)p  V  .  The  high  pressure  may  be  of  concern  if  the 
wall  can  only  withstand  a  limited  inpact  force. 


5.2  Qie-Fluid  Plow 

In  the  one-fluid  model  we  set  the  physical  properties  (density  and 
viscosity)  of  the  jet  fluid  equal  to  that  of  the  contaminant.  Then,  the 
interaction  between  the  jet  and  the  contaminant  becomes  a  flow  phenomenon 
that  takes  place  within  a  single  fluid.  To  justify  this  treatment,  we  only 
examine  the  viscosity  effect  because  the  densities  of  water  and  contaminant 
are  approximately  equal.  For  this  purpose  we  computed  several  flow  cases 
using  the  two-fluid  model  and  tested  the  sensitivity  of  the  results  to  the 
viscosity  of  the  jet  fluid  (water). 

Figure  26  presents  the  flow  patterns  resulting  from  two  viscosity 

values  of  jet  fluid,  v  *  0.98  mm^/s  (water)  and  v  *  98  mm^/s  (same  as  the 

w  w 

contaminant  viscosity).  Figure  27  shows  the  same  cases  when  the  water  layer 
over  the  droplet  is  increased  from  0.2  mm  (used  in  Figure  6)  to  1  mm.  The 
flow  patterns  in  the  left  and  the  right  hand  columns  are  very  similar  in  all 
cases.  Table  2  summarizes  the  displacements  S  and  the  mean  velocities  S  of 
droplet  upstream  edge  as  functions  of  time  for  the  flow  shown  in  Figure 
26.  Despite  a  drastic  variation  in  the  viscosity  of  the  jet  fluid,  the 
difference  in  the  mean  velocities  *S  is  only  of  the  order  of  10  -  15%.  A 
difference  of  this  magnitude  is  not  important  because  the  mean  velocity  at 
the  early  stage  (at  0.15  ms)  is  so  large  (S  >  2.7  m/s,  see  Table  2)  that 
practically  the  droplet  is  displaced  almost  instantly  after  application  of 
the  jet.  The  insensitivity  of  the  flow  to  the  viscosity  of  jet  fluid 
justifies  the  use  of  the  one-fluid  treatment.  (bntrary  to  the  jet 

viscosity,  the  contaminant  viscosity  is  important  as  presented  in  the 
following.  It  should  be  noted,  however,  that  this  flow  treatment  is  not 
suitable  if  the  impact  point  of  the  jet  is  located  far  away  from  the 
droplet.  In  that  case,  a  thick  viscous  layer  will  develop  in  front  of  the 
droplet  causing  a  far  shorter  displacement  of  the  droplet  than  by  a  less 
viscous  jet. 

5.2.1  Plow  Patterns  and  Viscosity  Effects.  Figure  28  presents  a 
sequence  of  flow  developments  following  the  initiation  of  the  jet 
impingement  in  Figure  lb.  The  stream  first  spreads  out  on  the  wall  and  then 
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engages  the  droplet  and  finally  is  lifted  off  the  wall.  The  profile  of  the 
jet -contaminant  interface  cannot  be  identified  in  the  figure  as  a  result  of 


Figure  24.  Mean  Velocity  of  Droplet  Upstream  Edge,  S,  vs.  Time 

after  (bmmencement  of  Jet  Flow,  t,  for  Various 

Combinations  of  Jet  Velocity,  Vj,  and  Jet  Size,  Dj 

the  one-fluid  treatment.  However,  by  embedding  marker  particles  in  the 
droplet  region  we  are  able  to  track  the  interface.  The  map  of  the  marker 
particles  in  the  case  of  Figure  28  is  presented  in  the  first  column  of 
Figure  29.  Other  columns  correspond  to  Increased  fluid  viscosities. 
Evidently,  the  viscosity  has  a  strong  effect  on  the  profile  and  the 

displacement  of  the  droplet.  The  viscosity  does  smoothen  out  the  interface 

and  retards  the  movement  of  the  contaminant.  Mother  interesting  phenomenon 

we  have  found  is  that  even  if  the  contaminant  viscosity  is  as  small  as  9.8 
mm^/s  (i.e.,  10  times  water  viscosity),  the  downstream  end  of  the  droplet 

remains  unchanged  in  shape  and  location  until  a  half  of  the  droplet  in  the 
front  has  been  broken  up  or  distorted. 
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Figure  25.  Pressure  Distribution  on  Impingement  Wall 
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Figure  26.  Flow  Patterns  (Two-Fluid  Flow,  Thickness  of  Water  layer 

above  Droplet  "0.2  mm,  vc  ■  98  mm^/s,  Vj  *  10  m/s,  0  *  45 


v  •  0.98  mm  Is 


[fttlttttttttllil 


t  ■  0. 003  ms 


I  nitia I  Location-*  t-  0  05  ms 


!!!>!!!!!! 


*98  mm  /s 

W 


HHIHU 


! ! ; ! : ! jj i : s: j 


rmsu 


t  ■  0. 10  ms 


1.1  II  mill.. 

::::::: 

iiiIiiiiiiiiiiSiiigiiiiiilllinll??l  j  j  I  ill! 


*  •  ♦  ♦ 
*'  ‘  •  ttSISSWWVNNW 
•••••Mini 
!!!?••"*  •  •  •  ••••*•!— 


•  ••••••••••  ••••  *  *•••••••  *  • 

Liiliiillliiliiiii  i  i  i  i  i  i  i  i  i  ?  i  I  *  iU*';  i i  I 


t  •  0. 15  ms 


Figure  27.  Flow  Patterns  (Two-Fluid  Flow,  Thickness  of  Watero 
Layer  above  Droplet  =  1  mm,  v  =  98  mm^/s,  6  =  45  ) 
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Figure  28.  Flow  Development  (One-Fluid  Flow,  V,  =  10  m/s, 
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Table  2.  Displacement  S  (mm)  and  Mean  Velocity  S  (m/s)  of  Droplet 
Upstream  Edge  vs.  Time  after  Cbmmencement  of  Jet  Flow  for 
Real  and  Ficticious  Viscosities  of  Jet  Fluid  v 
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Figure  29.  Evolution  of  Droplets  Corresponding  to  Various 
Viscosities  (One-Fluid  Flow,  6  *  45°) 

5.2.2  Effects  of  Jet  Velocity  and  Jet  Site  on  Jet  Performance .  Figure 
30  presents  the  displacements  S  versus  time  for  the  fluid  viscosity  vc  - 
v  -  98  mm^/s  and  the  angle  of  jet  incidence  0  ■  45°.  The  jet  stream  at  a 
higher  velocity  engages  the  droplet  earlier  and  thus  displaces  its  upstream 
edge  farther  at  a  given  time  after  the  impingement.  If  the  time  is  measured 
from  the  moment  at  which  the  displacement  starts,  then  the  origins  of  all 
three  curves  in  Figure  30  coincide  as  shown  in  Figures  31  and  32.  hi  these 
figures,  the  three  dashed  lines  are  derived  from  the  three  solid  lines  and 
represent  mean  velocities  S  of  the  droplet  upstream  edge  as  functions  of  S 
and  t.  They  show  that  an  increase  of  the  jet  velocity  greatly  Improves  the 
jet  performance. 

Figures  33  and  34  show  flow  patterns  for  8  -  0°  and  two  jet  sizes,  Dj  - 
D  and  Dj  *  2.5  D,  respectively,  where  D  **  0.6  mm,  same  as  the  height  of  the 
droplet.  We  observe  that  the  additional  amount  of  jet  fluid  from  a  larger 
jet  size  does  not  provide  a  significant  gain  in  droplet  displacement,  but 
simply  f lows  over  the  droplet.  We  have  not  conducted  investigations  of 
cases  with  9  >  0°,  but  it  can  be  reasonably  expected  that  at  the  same  jet 
flow  rate,  a  jet  with  a  higher  velocity  will  outperform  a  jet  with  a  larger 
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Figure  32.  Displacement  S  and  Mean  Velocity  S  of  Droplet  Upstream  Edge 
vs.  Time  after  Displacement  S  takes  place 


Figure  33.  Comparison  of  flow  Developments  between  Two  Jet  Sizes,  Dj 
=  D  and  2.5D  (One-Fluid  Flow,  v  =  V  =9.8  mm/s,  V 
=  10  m/s,  9  =  0°,  D  =  0.6  mm)  W 


size.  This  is  the  same  conclusion  as  in  the  case  that  the  droplet  was 
initially  covered  with  water. 

5.2.3  Pressure  Distribution  on  Impingement  Wall.  The  computed 
pressure  distributions  on  the  wall  for  various  times  are  shown  in  Figure 
35.  Their  maximum  value  is  twice  the  steady-state  stagnation  pressure  of 
the  jet.  This  is  much  lower  than  the  pressure  shown  in  Figure  25  in  which 
the  droplet  was  initially  covered  by  a  water  layer. 

VI.  JBT- CONTAMINANT  INTERACTIONS  IN  CONFINED  GEOMETRIES 

In  the  following  we  will  briefly  examine  flow  interactions  taking  place 
in  confined  geometries  other  than  the  plane  wall.  The  geometries  considered 
are  cavities  in  a  wall  and  corners  of  two  perpendicular  walls  shown  in 
Figures  36  and  37.  The  basic  difference  between  the  cavity  and  the  corner 
is  that  in  a  cavity  there  is  no  flow  through  the  left  boundary  of  the  flow 
region.  In  each  geometry,  two  values,  0  »  45°  and  0  90°,  were  used  for  the 
angle  of  jet  incidence.  The  contaminant  in  the  confined  geometries 
initially  can  be  in  form  of  a  droplet  or  covers  the  entire  bottom  surface 
with  a  thin  water  layer  above  it.  The  two-fluid  model  was  used  for  flow 
calculation.  We  now  present  the  computer-generated  flow  patterns  and 
typical  pressure  distributions  on  the  bottom  surfaces  of  the  geometries. 

6.1  Flow  Patterns. 

6.1.1  Interactions  In  Cavities.  Figure  38a  shows  the  flow  generated 
by  a  water  jet  impinging  on  the  upper  surface  of  the  cavity  filled  with 
water  and  a  contaminant  droplet  located  at  the  corner,  as  depicted  in  Figure 
37a.  All  of  the  following  conputer-generated  plots  have  been  magnified 
three  times  in  the  vertical  direction.  The  plots  in  the  left  and  the  right 
columns  correspond  to  the  45°  and  90°  impingement,  respectively.  The  figure 
shows  that  the  angle  of  Incidence  has  a  profound  effect  on  the  flow 
direction  and  accordingly,  the  movement  of  the  droplet.  In  the  45°- 
lmpingement,  the  main  stream  tends  to  move  toward  the  wall  on  the  other  side 
of  the  cavity  before  it  exits  the  region,  while  in  the  90°-impingement,  a 
large  fraction  of  the  jet  stream  exits  the  cavity  right  next  to  the  entrance 
region  of  the  jet  stream.  At  0.3  ms  there  is  a  small  amount  of  contaminant 
stagnating  at  the  right  corner  in  the  90°  case,  but  not  in  the  45°  case. 
The  jet  impinging  at  45°  seems  to  have  more  cleaning  power  for 
decontamination  of  the  cavity.  There  are  two  vortexes  seen  in  the  flow 
field,  one  adjacent  to  the  entrance  of  the  jet  and  the  other  inside  the 
droplet.  The  second  vortex  can  be  clearly  identified  in  Figure  38b. 
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Figure  34.  Evolution  of  Droplets  Corresponding  to  Figure  32 


Figure  35.  Pressure  Distributions  on  Impingement  Wall 
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Figure  38b.  Evolution  of  Droplets  Corresponding  to  Figure  38a 


If  a  droplet  Is  Initially  located  at  the  left  end  corner  of  the  cavity 
as  depicted  In  Figure  36b,  then  again  the  jet  in^inging  at  45°  provides  a 
higher  rate  of  cleaning,  as  shown  in  Figures  39a  and  39b. 

Figures  40a  and  40b  show  another  case  that  the  impingement  takes  place 
in  the  central  part  of  the  cavity.  In  the  90°-inpingement ,  the  main  stream 
of  the  jet  does  not  extend  far  enough  to  the  vertical  walls,  causing  only 
little  movement  of  the  contaminant  along  the  wall. 

Figure  41  shows  flow  developments  corresponding  to  the  configuration  in 
Figure  36d  in  which  the  cavity  is  initially  filled  with  contaminant  covered 
by  a  thin  water  layer.  Ihe  jet  impinging  at  90°  appears  to  have  more  power 
to  displace  the  contaminant  at  very  early  times.  However,  this  superiority 
disappears  after  t  >  0.18  ms  as  a  result  of  a  faster  development  of  viscous 
layer  along  the  vertical  wall  for  v  ■  90°. 
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Figure  39a.  Flow  Developments  Corresponding  to  Figure  36b 


6.1.2  Interactions  at  Corners  of  Two  Perpendicular  Walls 

Figures  42  and  43  show  flow  developments  corresponding  to  the 
configurations  of  Figures  37a  and  37b,  respectively.  Both  45°  and  90° 
result  in  similar  flow  patterns  except  that  there  is  a  small  amount  of 
contaminant  stagnating  at  the  corner  in  the  90°-impingeiaent. 

In  summary,  we  conclude  from  the  flow  patterns  presented  that  an 
Inclined  jet  impingement  is  more  effective  for  decontamination  in  various 
confined  geometries. 


6.2  Pressure  Distribution  on  Bo t toe  Surfaces  of  Confined  Geometries 


Figures  44  through  47  present  the  instantaneous  pressure  distributions 
on  the  bottom  walls  in  the  flow  cases  discussed  above.  In  any  flow  case, 
the  peak  pressure  far  exceeds  the  corresponding  steady-state  stagnation 
pressure  of  the  jet.  Uith  the  same  contaminant  configuration,  the  pressure 
rise  is  relatively  higher  in  cavities  than  at  corners  as  a  result  of  more 
flow  confinements  in  cavities.  In  general,  the  pressure  rise  in  a  confined 
geometry  is  a  function  of  the  dimension  ratio  of  the  jet  to  the  geometry, 
the  amount  of  contaminant  originally  resided,  the  contaminant  viscosity,  the 
location  of  jet  impingement ,  the  angle  of  jet  incidence,  and  so  forth. 


VII.  GENERAL  REMAKS 

The  following  general  remarks  can  be  made  based  on  the  results  reported 
in  the  previous  chapters. 


7.1  Application  of  Air  Jets 

A  great  advantage  of  decontamination  by  air  jets  is  that  there  is  no 
shortage  in  the  supply  of  jet  fluid  (air).  However,  air  jets  are  less 
effective  in  conparison  with  water  jets  because  of  the  much  lower  air 
density  which  produces  a  much  smaller  impingement  force  on  the 
contaminant.  Because  air  density  is  only  0.0012  of  water  density  in 
atmosphlcal  pressure,  the  air  speed  has  to  be  286  m/s  in  order  to  produce  a 
dynamic  pressure  equivalent  to  that  of  water  jets  at  10  m/s.  The  jet  at 
such  a  high  speed  may  constitute  a  noise  hazard.  In  addition,  often  there 
are  contaminant  stains  left  on  the  surface  cleaned  by  air  jets. 


7.2  Arrays  of  Jets 

In  Chapter  V  we  concluded  that  at  a  given  jet  flow  rate  a  small  high¬ 
speed  jet  provides  a  larger  cleaning  speed  and  consumes  less  jet  fluid  than 
a  large  low-speed  jet.  To  clean  a  large  area,  it  is  reasonable  to  assemble 
a  number  of  small  high-speed  jets  in  arrays  to  preserve  their  superior 
performance. 


WATER  LAYER 


There  are  many  possible  arrays  of  single  jets.  Depicted  in  Figure  48 
are  some  general  arrangements.  Single  jets  at  an  Inclined  angle  may  be 
assembled  in  form  of  straight  lines  "a,"  concaves  “a,"  a  circle  "c,  "  or 
convexes  “d."  The  arrangement  "a"  is  effective  when  the  array  is  moved 
forward,  but  not  sideways.  The  arrangement  "b”  will  push  the  contaminant  to 
the  central  region  and  has  to  carry  it  all  the  way  when  the  spray  is  moved 
forward.  It  may  perform  better  than  "a”  when  moved  laterally.  The 
arrangement  "c"  can  provide  effective  cleaning  in  all  directions.  The 
arrangment  “d”  may  have  a  performance  superior  to  “c"  because  the  impact 
locations  of  the  jets  in  the  two  rows  in  "d"  can  be  arranged  in  an  alternate 
pattern  and  hence  each  jet  in  the  array  can  be  used  more  efficiently. 
Experiments  are  needed  to  test  this  advantage. 


7.3  Justification  of  Validity  of  Two-Dimensional  Plow  Model 


The  flow  from  a  single  round  jet  impinging  on  a  surface  is  three- 

o 

dimensional.  If  the  angle  of  jet  incidence  is  90  and  the  surface  is  flat 
and  smooth,  the  fluid  uniformly  spreads  out  in  all  directions.  When  the  jet 
is  Inclined  gradually,  the  stream  in  the  lateral  direction  shrinks  and  the 
stream  in  the  front  grows  accordingly.  Eventually,  the  flow  pattern 
resembles  that  formed  by  a  two-dimensional  jet.  Also,  when  a  number  of 
round  jets  separated  by  a  small  distance  from  each  other  are  arranged  in  a 
straight  line  or  in  an  arc  pattern,  the  resultant  flow  on  a  flat  surface 
resembles  a  two-dimensional  flow. 
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Figure  48.  Arrays  of  Jets 


In  applications,  the  jet  spray  will  be  composed  of  a  number  of  single 
jets  assembled  as  shown  In  Figure  49  and  the  jet  flow  will  be  directed  at  an 
impingement  angle  less  than  90°.  Thus,  the  two-dimensional  models  developed 
in  Chapter  II  are  suitable  for  an  investigation  of  fundamental  flow 
characteristics  of  the  interaction  between  such  arrays  and  contaminant 
droplets. 


7.4  Rilsatlc 


A  jet  with  pulsating  flow  creates  a  great  transient.  Each  pulse  of  the 
flow  acts  as  an  initiation  of  a  new  jet  flow  and  the  pulsation  produces  a 
series  of  impact  pressures  on  the  surface.  These  impact  pressures  are  much 
greater  than  the  pressures  produced  by  a  continueous  jet  flow  without 
pulsation,  as  has  been  shown  in  Chapter  V.  A  large  ingtact  pressure  is 
particularly  helpful  to  break  up  and  displace  hard-to-remove  contaminant 
droplets.  The  flow  interaction  involved,  however,  is  extremely  complex  and 
it  is  not  clear  how  the  interaction  is  related  to  the  frequency  of  pulsation 
and  to  other  flow  parameters.  During  the  period  between  two  adjacent 
pressure  pulses,  there  is  a  little  or  no  force  acting  on  the  contaminant 
droplet.  It  is  possible  that  the  contaminant  ceases  to  move  or  even  moves 
backward  if  the  water  layer  which  covers  it  is  thick  and  if  the  frequency  of 
pulsation  is  not  appropriate. 


We  have  attempted  to  characterize  this  flow  using  the  SOLA-VOF  code, 
but  were  unable  to  obtain  reasonable  results.  The  difficulties  encountered 
were  the  numerical  instability  and  divergence.  A  modification  of  the 
numerical  method  employed  In  the  code  possibly  may  resolve  the  problem. 


Simulation  of  a  Movii 


In  this  report  we  have  assumed  that  the  jet  is  stationary  in  the 
physical  space.  Simulations  were  also  attempted  for  a  jet  moving  toward  the 
contaminant  droplet,  but  they  were  not  sucessful  because  of  numerical 
problems.  Mew  code  development  is  necessary  in  order  to  adapt  the  SOLA-VOF 
code  to  this  problem. 


VIII.  SUMMARY  AND  CONCLDSIONS 

TVo  flow  models,  called  one-fluid  flow  and  two-fluid  flow,  have  been 
used  to  investigate  jet-contaminant  interactions  on  a  plane  wall.  Both 
models  represent  two-dimensional  viscous  flows,  governed  by  the  Navier- 
Stokes  equations.  The  computations  were  done  with  the  SOLA-VOF  finite 
difference  code  in  which  we  implemented  a  viscosity  relation  to  adapt  the 
code  to  handle  the  present  case  that  the  two  fluids  involved  have  different 
viscosities. 


The  report  presents  computer-generated  flow  patterns  which  show  typical 
flow  developments  following  the  jet  Impingement.  It  is  found  that  the 
contaminant  viscosity  has  strong  effects  on  the  displacement  and  the  profile 
of  the  droplet.  In  the  high  viscosity  case,  a  thick  viscous  layer  develops 
quickly  inside  the  contaminant  region  and  the  jet  stream  is  lifted  off  the 
wall.  As  a  result,  the  contaminant  with  high  viscosity  is  much  more 
difficult  to  remove.  (bmputed  results  show  that  the  angle  of  jet  incidence, 
the  jet  velocity,  and  the  cross-sectional  area  of  the  jet  are  inportant 
parameters  that  govern  the  performance  of  the  jet.  The  results  also  show 
that  the  jet  can  perform  in  the  most  effective  and  most  efficient  way  in 
decontamination  at  an  angle  of  incidence  in  the  range  of  45®-60°.  Also,  it 
is  more  advantageous  to  use  a  jet  spray  composed  of  a  number  of  small  high¬ 
speed  jets  than  one  consisting  of  a  single  large  low-speed  jet. 

Jet  contaminant  interactions  taking  place  in  confined  geometries,  such 
as  cavities  and  corners  of  two  perpendicular  walls,  were  also  examined.  The 
results  show  that  an  inclined  jet  is  more  effective  than  a  normal  jet  in 
confined  geometries. 

In  all  flow  cases  studied,  the  Instantaneous  Impact  pressure  on  the 
wall  can  far  exceed  the  corresponding  steady-state  stagnation  pressure  of 
the  jet 
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area  which  has  been  decontaminated 


sound  speed,  Eq.  (1) 
a  constant,  I>0.6  mm 


thickness  of  jet  (jet  diameter  in  three-dimensional  cases) 
"fractional  volume  of  fluid"  function,  Eq.  (4) 


pressure 


pressure,  Eq.  (1) 
variable,  Eq.  (8) 
jet  flow  rate,  m^/s 

displacement  of  droplet  upstream  edge,  mm 
mean  velocity  of  droplet  upstream  edge,  m/s 


a  constant.  Figure  5 


U  fluid  velocity  on  the  wall 

u  fluid  velocity  component,  Eq.  (1) 

Uj  fluid  velocity,  Eq.  (7) 

v  fluid  velocity  component,  Eq.  (1) 


V  volume  of  jet  fluid  consumed 


Vj  jet  velocity 
x  coordinate,  Eq.  (1) 


x j  distance,  Eq.  (7) 
y  coordinate,  Eq.  (1) 


0  angle  of  jet  incidence 


k*  o • 


NGMENCLATUKE 


v  kinematic  viscosity  of  fluid  mixture,  Eq.  (1) 

v  kinematic  viscosity  of  fluid  mixture,  Eq.  (1) 

c 

v  kinematic  viscosity  of  water 
w 

p  density  of  fluid  mixture,  Eq.  (1) 

p  density  of  contaminant 

c 

Pw  density  of  water 

Subscripts 

c  contaminant 

j  jet 

w  water 
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